Cr02: a self-doped double exchange ferromagnet 
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Band structure calculations of Cr02 carried out in the LSDA+U approach reveal a clear picture 
of the physics behind the metallic ferromagnetic properties. Arguments are presented that the 
metallic ferromagnetic oxide Cr02 belongs to a class of materials in which magnetic ordering exists 
due to double exchange (in this respect Cr02 turns out to be similar to the CMR manganates). 
It is concluded that Cr02 has a small or even negative charge transfer gap which can result in 
self-doping. Certain experiments to check the proposed picture are suggested. 
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There has been a revival in interest in 3d transition 
metal oxides during the last decade. This was initially 
stimulated by the discovery of High-Tc superconductivity 
in complex copper oxides, and more recently by the ac- 
tive study of the colossal magnetoresistance manganates 
(CMR) Lai_^(Ca,Sr)a.Mn03. These latter systems in 
the most interesting composition range are metallic fer- 
romagnets. This is interesting in itself because ferromag- 
netic ordering is rare among the oxides: most of them 
are antiferromagnetic or ferrimagnetic with dominating 
antiferromagnetic interactions. We need to meet certain 
special conditions to stabilize ferromagnetism. One of 
the main mechanisms invoked to explain ferromagnetic 
ordering in these systems is the double exchange mecha- 
nism 0, although it is not the only one 

Another very well known ferromagnetic and metallic 
compound is chromium dioxide Cr02 widely used in mag- 
netic recording tapes. In its formal A+ valence state 
Cr has two 3c? electrons in t2g orbitals which in a sim- 
ple picture of strong correlations would suggest a Mott 
insulating-like ground state with S—1 local moments and 
most likely antiferromagnetic spin order. This seems to 
be about as far from the actual observed properties as 
one can get. In this letter we will address this problem 
using band structure methods supplemented with local 
Coulomb and exchange interactions (LSDA-I-U, [||) inter- 
preted in terms of local electronic configurations. We will 
argue that the d electrons can be divided into a localized 
"core" of spin ^ and an itinerant d electron propagating 
through these "cores" resulting in a double exchange- 
like mechanism for the ferromagnetic order much as in 
the manganates. We also show that strong electron cor- 
relation effects do not in this case lead to an insulating 
ground state for reasonable values of the d-d Coulomb in- 
teraction and perfect stoichiometry because Cr02 should 
be viewed as a small or even negative charge transfer gap 
M material in the Zaanen-Sawatzky- Allen (ZSA) scheme 



quite unlike the parent CMR material LaMnOa which 
is a Mott-Hubbard insulator. This leads for Cr02 quite 
naturally to a phenomenon which could be referred to as 
self-doping resulting in a non-integral 3d band occupa- 
tion §. 

Cr02 is a ferromagnetic metal with a saturation mag- 
netic moment of 2.00 fiB and a low temperature resistiv- 
ity with a nearly T^ temperature dependence ||^. Band 
structure calculations in LSDA |^ explain this behaviour 
as that of a half-metallic ferromagnet with a gap in the 
minority spin band resulting in the integral magnetic mo- 
ment per formula unit in spite of the strong covalency ef- 
fects. However the magnetic susceptibility in the param- 
agnetic phase exhibits a Curie- Weiss-like behaviour with 
a local moment of also 2 fj,B indicating that the mecha- 
nism for ferromagnetic behaviour is not a band (Stoner- 
like) mechanism. Local moments indicate strong correla- 
tion effects and this was in fact suggested by early pho- 
toemission data [|j which looked more like those of a 
semiconductor with a vanishing density of states (DOS) 
at the Fermi energy and a spin polarization of nearly 
100 % for binding energies about 2 eV below the Fermi 
level. To resolve this controversy we carried out the cal- 
culations in the LSDA+U approach which is superior to 
LSDA since it can indeed yield a gap if the local Coulomb 
interactions are large enough. 

Cr02 has a rutile structure (space group D^^: 
P42 / ninm) in which the unit cell consists of two formula 
units. The Bravais lattice is tetragonal (c/a=0. 65958) 
with a lattice constant of a—AA21 A The Cr atoms 
form a body-centered tetragonal lattice and are sur- 
rounded by distorted oxygen octahedra. The octahedra 
surrounding Cr at the body's center and corner positions 
differ by a 90° rotation about the c axis (see also Fig.|l| 
where Crl and Cr2 atoms are in the corners of the unit 
cell and Cr3 atom is in the body center position). 
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FIG. 1. Angular distribution of xy (left) and (yz + zx) 
(right) electron spin density for the nearest Cr neighbors. Di- 
rection of crystallographic c-axis is shown by arrow. Solid 
circles denote oxygen atoms. 

The simplest CrOg cluster calculations showed, that 
this kind of distortion of the octahedra (elongation along 
c axis) leads to the new natural basis for the t2g orbitals: 
xy, {yz + zx) and [yz — zx) in a local coordinate system 
(LCS) for every octahedron. In this LCS the z-axis is 
directed to the apex oxygen and the x and y axes - to 
the basal plane oxygens. Again refering to the Figj|, in 
which for Cr3 atom z-axis is directed to the 03 oxygen 
and X, y axes - to the 01 and 02 oxygens. At the same 
time, the 01 and 02 become apex oxygens for Crl and 
Cr2, and z axis is directed to these oxygens there. We 
will demonstrate below that this distorted structure is 
responsible for rather peculiar properties of the d bands. 

LSDA+U calculations were performed in the linearized 



muffin-tin orbitals (LMTO) approach |11| with atomic 
spheres radii of 2.06 a.u. for Cr and O and 1.78 and 
1.62 a.u. for the two kinds of empty spheres used (8 
per unit cell in total) to fill the empty space as much 
as possible. The spherical harmonics were expanded to 
the value of lmaa;=3 and 2 for the atomic and empty 
spheres, respectively. For the Brillouin zone (BZ) inte- 
gration 1300 k-points were used. The screened U and J 
parameters used in the LSDA-I-U scheme were calculated 
by constrain method taking into account participation 



of eg electrons in screening of the t2g electrons [|2j and 
found to be 3 and 0.87 eV, respectively. 

The results of this calculation are shown in Figs. ^jj. 
First, for U=3 eV the material is found to be a half- 
metallic ferromagnet, similar to the previous LSDA re- 
sults ||]. However, inclusion of the electron correlations 
modifies the electronic structure. Together with the shift 
of the minority spin DOS from the bottom of the con- 
duction band to higher energy, there appears a dip in 
the majority spin DOS at the Fermi level which can be 
taken as an indication of the tendency toward a gap for- 
mation. The decrease of the DOS at E^? in comparison 
with that obtained in Q may be partially responsible for 
the reduced signal close toEp observed in photo-electron 
spectroscopy In our calculations, however, Cr02 is 
still a metal albeit with the reduced DOS at the Fermi 
level. An explanation of this discrepancy may lie, e.g., 
in the modified situation at the surface. Note also that 
more recent photoemission data show a finite DOS at 
Fip, and the spectrum obtained in [p^ is in rather good 
agreement with our calculations. 

The tendency to open a gap in an energy spectrum 
which we noticed above was confirmed by calculation 
with larger values of U. It was found that the gap would 
indeed open at U>6 eV which would make Cr02 a Mott- 
Hubbard insulator. 

There are other aspects of the calculation that are very 
interesting and may provide the clues to the basic mecha- 
nism for the ferromagnetic metal behaviour. First notice 
the almost dispersionless majority spin band at about 
1 eV below Ej? in Fig.^ over a large region of the BZ 
which as is indicated by the black circles has almost pure 
d character. This is observed as a peak in the DOS in 
Figs.^ and ^, and corresponds to strongly localized xy 
orbitals completely occupied by one majority spin elec- 
tron. Because these states have almost pure d character 
they are very sensitive to the size of U and shift further 
away from Ej? as U increases. 




FIG. 2. LSDA-f U band structure of Cr02 for majority spin-sublattice in the vicinity of the Fermi level. Open circles denote 
the contribution of d states less than 25 %, light gray circles ~ between 25 and 50 %, gray - between 50 and 75 %, and black - 
predominantly d character of the band. Fermi level is denoted by horizontal solid line. 
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3. Total (per formula unit and both spins, a) and 
Cr 3rf- (b) and O 2p- (c) density of states of Cr02. 



The other t2g d-bands are strongly hybridized with O 
and form dispersive bands crossing the Fermi level. One 
can see from Fig.^ that predominantly the (yz + zx) is 
occupied at each site, see also Fig.|l|. An admixture of 
oxygen states strongly softens the influence of U on these 
bands and that is the reason that such a high value of 
U=6 eV is required for a gap to open in these bands. This 
indicates that we are dealing with a system of localized 
S=| xy electrons exchange coupled to broader band-like 
and strongly hybridized {yz + zx) states. In our case the 
bandwidth of d-p 7r-bands is comparable with the Hund's 
rule intraatomic exchange (^1 eV). This is reminiscent 
of a double exchange ferromagnetic mechanism of Zener 
101 . This is also consistent with the half- metallic nature 
of this material. Thus, the (half-)metallic character and 
ferromagnetism are indeed closely connected. 

An extra confirmation of this conclusion came from the 
calculation, in which we artificially imposed an antifer- 
romagnetic ordering. The simplest two-sublattice anti- 
ferromagnetic structure was chosen, with the spin in the 
body center sites opposite to the spins at the corners of 
the unit cell. The band structure obtained has a small 
gap ~0.15 eV exactly at the Fermi-surface. That means 
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FIG. 4. Partial Cr t2g majority spin density of states in the 
vicinity of the Fermi level. 

that Cr02 in an antiferromagnetic phase would have been 
an insulator. Thus, the metallicity and ferromagnetism 
in Cr02 do indeed support each other, which is in full 
agreement with the double exchange picture. 

The reason for this dualistic behaviour of the t2g states 
i.e some localized and others dispersive, can be found in 
the distortion from an ideal rutile-type crystal structure 
[0. Our calculations show that in the ideal structure 
with the same length of all 12 0-0 bonds in a perfect oxy- 
gen octahedron, all the t2g bands would have the same 
width. In the real distorted structure the edge sharing 
nearest Cr neighbors (Crl and Cr2 in Fig.|l|) move apart 
decreasing the xy bandwidth and shifting its center of 
gravity to low energies in comparison with the centers 
of gravity of the other two t2g bands. And now with 
the inclusion of Coulomb interactions, this xy band be- 
comes the first candidate to be fully occupied. Just that 
happened in our LSDA-I-U calculation. This picture has 
much in common with the band structure, proposed by 
Goodenough |jl^ on phenomenological grounds. 

It is possible to visualize the d orbitals mentioned 
above (in particular, occupied orbitals) by plotting the 
angular distribution of d electron spin density for these 
orbitals (Fig.|l|): 
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where Qm,m' = n'^^rn' ~ '^m,m' '^^ d spin occupation 
number for the case of xy orbital or 2*2 matrix for the 
case of (yz + zx) orbital obtained in our self-consistent 
calculation and Ym{0, cj)) are corresponding spherical har- 
monics. Note that {yz + zx) orbitals of Cr atoms hy- 
bridize with each other via 7r-overlap with 01 and 02 
oxygens (right panel of Figj^). 

Yet another interesting feature is that there are bands 
of predominantly O 2p character crossing the Fermi level 
close to the F-point (Figj2|). The absence of hybridiza- 
tion of these with the d states indicates that they have 
symmetries such that this mixing is not allowed. The 
fact that they extend to the Fermi level is important. 
This means that these states can be used as electron or 
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hole reservoirs causing a non-integral occupation of the 
d bands. This may be called self-doping. This again is a 
reason for the metallic behaviour in spite of quite large 
U values. 

The situation when the almost pure O 2p bands cross 
the Fermi level indicates that such materials should be 
considered as lying in the small or even negative charge 
transfer gap region of the ZSA diagram. The negative 
charge transfer gap behaviour is not unexpected for ox- 
ides with transition metals in very high oxidation states 
because of the increased electron affinity of these ions. 
This conclusion can in principle be checked by the ls-2p 
x-ray absorption on oxygen; simultaneously certain spin 
polarization should be transferred to oxygen which would 
be very interesting to study experimentally. In our calcu- 
lations the oxygen polarization turned out to be opposite 
to the polarization of d-ions and equal to -0.09 /is per 
oxygen. 

One more conclusion of the picture obtained is that, as 
is common for the double exchange systems, one should 
expect a negative magnetoresistance in the paramagnetic 
region. Previous results have shown that the magnetore- 
sistance of Cr02 is rather small but these mea- 
surements were done only below room temperature. It 
would be interesting to extend these measurements to the 
region of T>Tc=392 K. As there is no extra "complica- 
tions" in Cr02 as compared to CMR manganates (no 
strong Jahn- Teller effects, no disorder, no formally dif- 
ferent valence ions), this study would reveal a pure dou- 
ble exchange contribution to magnetoresistance which, 
when compared with similar data for manganates, could 
help to discriminate between the contribution of different 
mechanisms in the latter. 

Summarizing, we carried out an extensive band struc- 
ture calculations of Cr02 in the LSDA-I-U approach. We 
confirmed the half-metallic nature of Cr02, albeit with 
the dip of DOS at the Fermi level. It is demonstrated that 
there exist in Cr02 two groups of c?-electrons with signif- 
icantly different properties: one of the two d-electrons 
of Cr'^^ is essentially localized, and another is 7r-bonded 
with the oxygen p-orbitals and forms a partially filled 
narrow band. The width of this band is of the same or- 
der as the Hund's rule intraatomic exchange interaction. 
The resulting picture may be interpreted as an indica- 
tion that the ferromagnetism in Cr02 is due to double 
exchange mechanism. An extra confirmation of this con- 
clusion comes from the fact that in an antiferromagnetic 
phase Cr02 would have been an insulator. The strong 
contribution of the oxygen p-states at the Fermi surface 



explains why this material is not a Mott insulator in spite 
of the large U values and the relatively narrow d-bands 
and also shows that Cr02 belongs to the class of com- 
pounds with small or even negative charge-transfer gap 
leading to a self-doping. Several experiments are sug- 
gested which could check the proposed picture and could 
help to isolate the double exchange contribution to mag- 
netoresistance. 
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